Dynamics of four- versus two-terminal transport through chaotic quantum cavities 
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We consider multi-terminal mesoscopic transport through a well-conducting chaotic quantum 
cavity using random matrix theory. Four-probe resistance vanishes on the average and is not affected 
by weak localization. Its fluctuations are given by a single expression valid for arbitrary temperature, 
ac frequency, non-ideal coupling of the contacts, and in the presence of floating probes; surprisingly, 
they are governed by the dwell time only. In contrast, the two-probe transport additionally depends 
on the RC-time, which is interpreted as a property of the measurement scheme. We also predict a 
universal mesoscopic distribution of the phase of transmitted voltage in an ac experiment. 
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Introduction. Mesoscopic transport and its sample- 
to-sample fluctuations are fundamentally important for 
understanding quantum effects in electronic transfer 
through nano-structures. However, when a voltage is 
measured between the source and drain in a two-terminal 
geometry, a big classical contribution to transport shad- 
ows quantum effects of order e'^/h [1, 2]. The visibility of 
interference effects can be enhanced by diminishing the 
width of the contacts to reservoirs, but this also leads to 
Coulomb blockade of electrons [2] . Another approach is 
to keep the contacts well-conducting and apply a multi- 
terminal scheme with a voltage through the sample mea- 
sured by two additional contacts, see Fig. 1(a). The thus 
defined four-terminal resistance R^t, well-known to ex- 
perimentalists, has some benefits over the two-terminal 
one, i?2t, which is more popular among theoreticians. 

The central advantage of i?4t is to exclude the contact 
properties. Indeed, classical resistances are local, and a 
four-terminal scheme removes the effect of the contacts 
and therefore probes the sample's properties only [3]. 
Quantum mechanically resistance is not local, and this 
scheme does reduce the role of contacts but they still af- 
fect the measurement. Nevertheless even for a coherent 
sample there are of course profound differences between 
i?2t and i?4t. The former determines Joule heating and is 
therefore always positive and symmetric to the magnetic 
field inversion. In contrast, the latter need not be pos- 
itive [4, 5] and experiments demonstrate its fluctuations 
around zero [6-8] and magnetic field asymmetry [9]. 

Non-universal fluctuations in i?4t were widely investi- 
gated in metallic diffusive samples, where small quantum 
effects on top of a large classical average depend on the 
probe locations [10-15]. On the other hand, in a chaotic 
cavity [16, 17] the exact positions of probes are irrele- 
vant, the classic voltage drop vanishes, and an experi- 
ment would measure quantum effects directly. Remark- 
ably, i?4t in this generic geometry remained unexplored 
except for a recent initial experiment [18]. 

What could one expect from R^t in a quantum dot? Its 



dc two-terminal conductance shows a weak-localization 
correction (WL) ^ /h due to an enhanced return prob- 
ability for an electron, and universal conductance fluctu- 
ations (UCF) ~ (e^//i)^; a magnetic field destroys WL 
and reduces UCF by a factor 2 [16]. These results for 
ballistic quantum point contacts (QPCs) can be gener- 
alized to include either (i) a dephasing (floating) probe 
coupling [19], or (ii) non-ideal QPCs [20], or (iii) a low- 
frequency ac setup [21]. In (iii) the averaged conduc- 
tance [21] and shot- noise [22] in their leading order de- 
pend only on the charge relaxation time (RC-time) trCj 
sensitive to the Coulomb interaction energy e^/2C of the 
capacitor. Its appearance is natural: at high frequen- 
cies Lo > I/trc the capacitor conducts better than the 
contacts. However, WL and UCF, as well as the third 
current cumulant, additionally depend on the dwell time 
Td an electron typically spends in the dot [21, 23, 24]. 
As a result, even in linear transport there is no unique 
time-scale for dispersion of various quantities. 

At first one expects R^t to share the features of two- 
terminal WL and UCF. However, these expectations are 
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FIG. 1: (a) A multi-terminal chaotic cavity with a gate (ca- 
pacitance C). A current through source-drain contacts 1,2 
and a voltage drop between 3,4 (or 1,2) determine the four- 
probe Rut ~ V34/I (or two-probe R2t = V12/I) resistances. 
Additional floating probes draw no current, (b) The meso- 
scopic distribution P{(j}) (rescaled by xuuTd) of phase (j) of 
the transmitted ac voltage as a function of frequency, for 
cjTd — 0.1, 0.3, 1, full, long- and short-dashed, respectively. 
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wrong: i?4t docs not have weak localization, and its fluc- 
tuations are insensitive to magnetic field. The difi'erences 
are even greater when the dynamics of resistance fluctu- 
ations is compared: imlikc i?2t governed both by and 
trC) R4t turns out to depend only on tj. Surprisingly, 
fluctuations of i?4t are much simpler than UCF, and a sin- 
gle formula covers the regimes (i-iii) for arbitrary floating 
probes, imperfections in contacts, and frequency ui. 

Universal resistance fluctuations. This Letter com- 
pares statistics of two- and four-terminal resistances in 
quantum dots with non-ideal coupling to reservoirs and 
probes, biased by ac voltages, see Fig. 1(a). To be spe- 
cific let the current / flow only through contacts 1,2 and 
measure the voltage drop between 3,4 (or 1,2) to find 
R4t = V34/I (or R2t = V\2/I)- Each contact j with Nj 
orbital channels is specified by a diagonal matrix V j of 
its channel transmissions, and we assume good dimen- 
sionless conductances of the probes 1-4, trFj 1. At 
arbitrary temperature T and measurement frequency lo, 
in first two orders in the total QPC conductance trF, the 
ensemble averaged resistances are given by 
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Vanishing (-R4t) is not sensitive to the presence of time- 
reversal symmetry (TRS), fi = 1, or its breaking by mag- 
netic field, j3 = 2. In contrast, (-R2t) is corrected by WL 
depending only on ra- Similarly, TRS does (not) affect 
fluctuations of R2\,{Ri\) around their averages: 
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Fluctuations of -R4t, Eq. (3), the main result of this paper, 
are valid for any couplings of floating probes, and imper- 
fections in the contacts, while Eq. (4) is derived only for 
ballistic QPCs in a multi-terminal dot. 

In an eff'ectively zero-dimensional dot universal Eqs. 
(1--4) do not depend on disorder and contact positions. 
It is interesting that for conductance there are no re- 
sults of comparable generality for WL and UCF . We 
conclude that both weak-localization corrections to the 
two-terminal resistance, Eq. (2), and fluctuations of the 
four-terminal one, Eq. (3), are completely insensitive at 
all frequencies to the RC-time and thus to Coulomb inter- 
actions. We next explain how these results are obtained, 
and discuss the dc fluctuations of resistances. Later we 
generalize our results onto an ac setup and discuss the 
role of the chosen measurement scheme in fluctuations. 



We calculate the universal distribution of the phase of 
transmitted voltage in a four-terminal ac experiment. 

The calculation. We consider a chaotic quantum dot in 
a multi-probe setup with M > 4 contacts and a gate with 
capacitance C, see Fig. 1 (a). Each contact j = 1, M 
leading to a reservoir with a voltage Vjuj at the fre- 
quency uj is characterized by a set of transmission values 
Tj of its (imperfect) channels. This diffusive/ballistic 
dot is in the universal regime, when direct trajectories 
arc absent, and chaos validates random matrix theory 
(RMT) [16, 17]. The Thouless energy of a closed dot is 
much larger than the temperature T, mean level spac- 
ing A = 2TTh'^/{m Area), the excitation energy hoj in the 
ac experiment, and escape rate trF • A/2tt. The total 
coupling to environment is good, trF 3> 1, and it speci- 
fies the electronic dwell-time Td = /i/(trFA) in the dot. 
The leading Coulomb interaction effect is accounted for 
by a uniform self-consistent potential of the dot, found 
from charge conservation [21, 25]. The terms leading to 
Coulomb blockade are unimportant due to 1/trF <C 1. 
We use RMT and diagrammatic technique with a small 
parameter 1/trF ^ 1 for an energy-dependent scattering 
matrix S{e) of a dot with imperfect contacts [20, 26, 27] 
to find transport statistics. Further we take h = e = 1 
and consider spin-degenerate electrons, Vs = 2, the role 
of spin-orbit in i?4t being discussed elsewhere [28] . 

Linear transport coefficients are found from the 
screened frequency-dependent degenerate (M+l) x {M+ 
1) conductance matrix g'^j ^ — dliuj/dVjuj, where 

9ij,u, = fe> - Y^kl 9il,uj9kj,uj/ {Y.kl 9kl,cj - iujC), 

=2.../..tr[l.l,-l,5t(,),,^(, + .^)] 

X (/(£) - /(£ + a;/27r))/a;, i, j, k,l = l, M,(6) 

using summation over k or I, X^^j^^o 9ki ui — ^ 25]. 
All probes except the current source 1 and sink 2, / = 
Ji = —I2, are set to voltages such that they draw no 
current. We shift all voltages by —V4U, and eliminate the 
4th row and column from g'^ , invert the rest and obtain 
the resistance matrix R. It is used to explore statistics 

of i?4t = i?31 — i?32 and i?2t = Ru — R12 — R2I + -^22- 

DC setup, oj = 0. First, let us consider dc transport 
and ballistic QPCs at T = 0, when any trF" = Nj in 
Eqs. (2,3) and the integral, Eq. (5), is reduced to 1. In 
this limit the measiircmcnts in a two-terminal dot arc not 
affected by the Coulomb interaction, and both (-R2t) and 
its fluctuations are diminished by the TRS-breaking [29]. 
In a multi-terminal setup, however, chaotic scattering re- 
sults in a random voltage drop between any pair of volt- 
age probes. As a consequence, the average R^ vanishes, 
(i?4t) = 0, and does not have a WL contribution. In- 
deed, i?2t = (ff33C/44 - 534543) • det R is symmetric with 
respect to 1 -f^ 2, and its WL is the first correction due 
to self-intersections of an electronic trajectory from 1 to 
2. On the other hand, i?4t = (531542 — 532541) • det R [4] 
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is anti-symmetric, and any loop contributes equally to 
both terms. Consequently, WL corrections are canceled, 
and they appear only if (-R4t) 7^ 0, e.g. in a quasi-ld dis- 
ordered structure [30] with a special probe arrangement. 

Similarly, the TRS-breaking does not affect the fluc- 
tuations of R4t as shown by Eq. (3), in contrast with 
Eq. (4) for R2t- Indeed, inverting time we must invert 
magnetic flux $ and swap the current- with voltage- 
probes. The ensuing symmetry i?2t(^) = i?2t(— gives 
9i?2t/9$ = at $ = 0, but there is no such a restric- 
tion on jR4t($). Fluctuations of i?2t(^') are analogous 
to a perturbation on a string with a free end, where the 
perturbation maximum corresponds to the maximum in 
i?2t-fluctuations at $ = 0. On the other hand, i?4t($) 
is analogous to an infinite string, so that zero field is as 
good as any other: neither TRS-breaking effects like WL 
in Eq. (2) or 2//3-coefficient of Eq. (4), nor "quenching" 
(diminished sensitivity to $ at $ = 0) appear for in 
chaotic dots [28]. In addition, in contrast to the classi- 
cal four-terminal scheme [3] , the properties of source and 
drain are never fully eliminated: a measurement affects 
the sample and so the mcsoscopic fluctuations in i?4t do 
depend on contacts. Importantly, in the experiment [18] 
the voltage probes are invasive. 

For the non-ideal contact coupling a tunneling limit 
for the voltage probes is often taken, trF^ <C 1; in the 
opposite limit, trFj, » 1, UCF presents a complicated 
6-th order polynomial of Tj [20]. In contrast, for a fixed 
set {Fj} the fluctuations of R^t given by Eq. (3) can be 
expressed only in terms of conductance trFj and Fano 
factor Tj, the ratio of shot noise to current [31], of each 
QPC. The result of a separate averaging over mesoscopic 
contact j, trF|/tr^Fj- (1 — (J'j))/(tr F^), is propor- 
tional to the QPC resistance with a numerical prefactor: 
it ranges from 1 for ballistic to 2/3 for a short diffusive 
wire, if disorder is increased; for a wide tunneling QPC 
it is equal to 1 if Tj <^ 1 arc the same, and 1/2 for a 
dirty interface [32] . We observe that, generally, a poorer 
coupling increases fluctuations of i?4t. 

Often conductors have more contacts then needed for 
a four-terminal measurement. A floating (unused) con- 
tact / does not to draw any current from the sample, 
but leads to decoherence [33] , its coupling is related to 
the inelastic scattering time Tin = h/AtrTf [19]. Data 
are often interpreted as an addition of several flctitious 
dephasing channels to the real contacts. However, here 
it is important to distinguish between a floating probe 
and a contact used for transport measurements. Simi- 
larly to Rcf. 19, we find that a fioating probe decreases 
fluctuations of i?4t by increasing the total coupling trF 
of the dot to reservoirs. Due to the zero-dimensional ge- 
ometry of the dot, they are independent of the number 
of such probes, and only the total trF/ is relevant. The 
voltage probes used to define i?4t and floating probes 
arc not equivalent, an increase of the former decreases 
fluctuations stronger. Indeed, in two measurements with 



unequal probes used as voltage and fioating ones and vice 
versa, the better the voltage probe conducts the smaller 
the fluctuations of i?4t. Importantly, if the dephasing 
probe model is used, dephasing channels should be in- 
cluded not into the real, but into the floating probes. 

AC setup, CO ^ 0. Now we generalize our dc re- 
sults to finite frequency lu to consider dynamics of resis- 
tances. The screened conductance in Eq. (6) becomes 
complex and depends not only on ra [34] , but also on 
Tjic = C/i/trF. Here is the electro-chemical capaci- 
tance, which relates the charge of the dot to the chem- 
ical potentials, it is defined by the Coulomb interaction 
strength, 1/Q, = l/C+A/i/^ [25]. A frequency 1/tr,c, 
when the capacitor and the contacts conduct similarly, 
is then a natural characteristic of a two-terminal mea- 
surement. The finite- frequency correlations and WL of 
conductance were considered only for ballistic QPCs at 
T,uj 0. Quantum effects in ac transport depended 
only on cjTd at low frequencies [21], but at higher uj they 
start to depend on wtrc as well. In contrast, for i?2t our 
Eq. (2) shows that at any frequency the WL is indepen- 
dent of trc; this is compatible with Tpc-dependent WL 
in conductance [21], because trc is eliminated by the 
matrix inversion. The reason is that WL only increases 
the return probability, but does not redistribute charge. 

Low-frequency conductance correlations can be gen- 
eralized to arbitrary frequencies u,u' using parameter 
z = — icjTRc/(l ~ i'-^Tji.c) and z' introduced similarly, 
and considering ballistic QPCs with rij = Ni/N for sim- 
plicity. For /3 = 2 we find the correlations, 

(glj.^gg;,^') ^ /•^^ (27r^T)2e-"/"^(e'"" - 1)(1 - e^"'") 
vlniUjUkni J rdww' sinh^(27r2 Tr) 

\ {kk/nk - 1 + zz'){5ji/ni -1 + zz') _ tz'^z'^ I 
^[ l-iTd{w + iv')/2 T^WW'J 

and for = 1 the same r.h.s. with k <H> I should be 
added. Results of Ref. [21] are reproduced by differen- 
tiation of Eq. (7) at co, co' = 0. Using Eq. (7) we 
find that fluctuations of i?2t in a dot with ballistic probes, 
N = Ni + N2+Nf, given by Eq. (4) manifestly depend on 
both Td and trc- Only the four-probe scheme can show 
if wtrc is an intrinsic scale of the dot itself, or if it comes 
from the chosen two-terminal measurement scheme. In 
other words, does the frequency-scale 1/trc survive in 
resistance statistics if the role of contacts is minimized? 
Equation (3) for Rit does not depend on tkc, and we 
conclude that this scale is extrinsic and can be under- 
stood as a contact property. Since Gaussian i?4t is fully 
characterized by Eqs. (1,3), we conclude that the quan- 
tum transport at arbitrary frequency oj is insensitive to 
screening only in the four-terminal scheme. 

Experimental relevance. Recently, Lerescu et al. [18] 
measured statistics of i?4t in dc-biased ballistic dots and 
compared with numerical data for Nf = 0. The rms R^t 
qualitatively agreed with the numerics, but was about 20 
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times smaller than they expected. These data at ^ 1 
indeed correspond to our predictions that (i?4t) = and 
its fluctuations are strongly diminished by an increased 
voltage probe coupling. Their numerical data are well- 
fitted by Eq. (3), which demonstrates a weak dependence 
on A^3 — A'4 and does not depend on /3 = 1,2. Reasons 
for the much reduced experimental arc strong dcco- 

herence, which might depend on couplings, and possibly 
direct source-drain trajectories, which result in voltage 
drop fluctuations smaller than expected by RMT. 

High-frequency experiments are more difficult. If a 
leakage current through C bypasses the drain, we have 
Ii^ 7^ —l2uj- Resistance i?4t relates the voltage drop V34 
to some linear combination of these currents, and our 
corresponds to — /2a;)/2. If the voltage drop is mea- 
sured via a transmission line with length comparable to 
27rc/w, the circuit should be also taken into account. The 
phase shift of transmitted voltage with respect to current 
source has a classical component due to the circuit and a 
fiuctuating part </> due the dot, and we consider its meso- 
scopic distribution If the line impedance between 

the exits 3,4 (1,2) is large compared to i?4t(-R2t), one 
has (j) = argi?4t. Using Eqs. (3,5) and (|i?4tp) calculated 
similarly, and normalizing § d(l)P{(f)) = 2'jt, we find 



and discussion. MB is supported by the Swiss NSF and 
MaNEP. 



1 



Lo -I^c [cxp(-2i(/.) J^J a + 



sm (0 - 0o) 



This TT-periodic distribution is insensitive to the probe 
properties and dephasing. The position (j)Q and the height 
of the distribution maximum arc defined by Eq. (8), and 
Fig. 1(b) presents P(0) for low T <C h/T^. The maxi- 
mum \/2/wTd at (po = WTd for wra ^ 1 is sharp and the 
reactive part of i?4t is small. At wtj ^ 1 the peak at 
00 = 7i'/2 — log2 WTd/(4a;Td) is small due to weak corre- 
lations between (re) active parts, and any phase ^ is pos- 
sible. For high T > fi uj,H/Td w hen RMT is still valid, 
one has P{(j)) = WTd/[-\/l -|- w^r^ — cos(2(/) — arctanwTd)], 
qualitatively similar to the low-T result. We expect that 
the universal distribution given by Eq. (8) can be found 
from phase measurements by varying either w or T. 

Conclusions. We discuss statistics of multi-terminal 
transport through a well-conducting chaotic quantum dot 
with imperfect coupling. Using RMT we find that four- 
probe resistance is unaff'ected by weak localization and its 
fluctuations arc governed only by the dwell-timc. Unlike 
two- probe transport, the four-terminal one is insensitive 
to the Coulomb interactions at any frequency, suggesting 
that the charge-relaxation time is connected to the mea- 
surement scheme. The fluctuations are given by a single 
analytical expression for arbitrary temperature, ac fre- 
quency, non-ideal coiipling of the contacts, and floating 
probes. We propose a universal mesoscopic distribution 
for the phase of voltage transmitted through the dot. 

Acknowledgments. Wc would like to thank K. Flens- 
berg, J. Gabelli, and C. M. Marcus for useful comments 



[1] 
[2] 

[s: 

[4] 
[5] 

[6: 

[8 

[9 
[10: 

[11 
[12 

[13 
[14 

[15 

[16 
[17; 

[18 
[19 
[20; 

[21 

[22; 

[23; 

[24 
[26 

[26 
[27 

[28 
[29 

[30 
[31 

[32 

[33 
[34 



Y. Imry, Introduction to Mesoscopic Physics (Oxford 
University Press, USA, 2002). 

B. L. Altshuler, P. Lee, and R. Webb, eds., Meso- 
scopic Phenomena in Solids (North-Holland, Amster- 
dam, 1991). 

L. J. van dor Pauw, Philips Rosorarch Reports 13, 1 
(1958). 

M. Biittiker, Phys. Rev. Lett. 57, 1761 (1986). 
M.Biittiker, IBM Journal Research Development 32, 317 
(1988). 

R. De Picciotto et al, Nature 411, 51 (2001). 

B. Gao et al, Phys. Rev. Lett. 95, 196802 (2005). 

G. Gunnarsson, J. Trbovic, and C. Schonenberger, Phys- 
ical Review B 77, 201405 (2008). 

A. D. Benoit et al., Phys. Rev. Lett. 57, 1765 (1986). 
M. Ma and P. A. Lee, Phys. Rev. B 35, 1448 (1987). 

B. Z. Spivak and A. Yu. Zyuzin, in [2]. 

H. U. Barangcr, A. D. Stone, and D. P. DiVincenzo, 
Phys. Rev. B 37, 6521 (1988). 

C. L. Kane, P. A. Lee, and D. P. DiVincenzo, Phys. Rev. 
B 38, 2995 (1988). 

S. lida and A. Miiller-Groeling, Phys. Rev. B 44, 8097 

(1991). 

B. Schmidt and A. Miiller-GroeUng, Phys. Rev. B 47, 

12732 (1993). 

C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

I. L. Aleiner, P. W. Brouwer, and L. I. Glazman, Physics 
Reports 358, 309 (2002). 

A. I. Lerescu, E. J. Koop, C. H. van der Wal, B. J. van 
Wees, and J. H. Bardarson, arXiv:0705.3179 (2007). 
P. W. Brouwer and C. W. J. Beenakker, Physical Review 
B 55, 4695 (1997). 

P. W. Brouwer and C. W. J. Beenakker, Journal of Math- 
ematical Physics 37, 4904 (1996). 

P. Brouwer and M. Biittiker, Europhysics Letters 37, 441 
(1997). 

M. L. Polianski, P. Samuclsson, and M. Biittiker, Physi- 
cal Review B 72, 161302 (2005). 

K. E. Nagaev, S. Pilgram, and M. Biittiker, Phys. Rev. 
Lett. 92, 176804 (2004). 

C. Petitjean et al, Phys. Rev. B 80, 115310 (2009). 
M. Biittiker, A. Prctre, and H. Thomas, Phys. Rev. Lett. 
70, 4114 (1993). 

P. W. Brouwer, Phys. Rev. B 51, 16878 (1995). 

M. L. Pohanski and P. W. Brouwer, Journal of Physics 

A: Mathematical and General 36, 3215 (2003). 

M. L. Polianski, Physical Review B 80, 241301 (2009). 

P. W. Brouwer, A. Lamacraft, and K. Flensberg, Physical 

Review B 72, 075316 (2005). 

T. Ihn et al, Physica E 40, 1851 (2008). 

Y. M. Planter and M. Biittiker, Physics Reports 336, 1 

(2000). 

K. M. Schep and G. E. W. Bauer, Phys. Rev. Lett. 78, 

3015 (1997). 

M. Biittiker, Phys. Rev. B 33, 3020 (1986). 

J. B. Pieper and J. C. Price, Phys. Rev. B 49, 17059 

(1994). 



